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Towards  Development  of  an  Improved  Technique  for  Remote  Retrieval 

of  Water  Quality  Components: 

An  Approach  Based  on  the  Gordon’s  Parameter  Spectral  Ratio 

Sokoletsky,  Leonid1;  Gallegos,  Sonia 

'Elizabeth  City  State  University  1704  Weeksvillc  Road,  Elizabeth  City,  NC,  27909,  United  States; 
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INTRODUCTION 

Remote-sensing  retrieval  of  water  quality  components  (WQCs)  such  as  ehlorophyll-a,  suspended 
solids  concentration,  or  dissolved  organic  matter  absorption  is  based  on  relationships  between  the 
WQCs  and  underwater  spectral  radiances  or  reflectances.  The  use  of  spectral  radiance  ratios  or  spectral 
reflectance  ratios  largely  decreases  errors  caused  by  the  angular  scattering  features  of  water  particles, 
geometrical  conditions  of  illumination  and  observation,  and  bottom  effeet.  Additionally,  it  flattens  the 
non-linear  behavior  of  reflectance  as  a  function  of  Gordon’s  parameter  G  =  bb/(a+bt>)  ( a  and  bb  are  the 
absorption  and  backscattcring  coefficients,  respectively).  Numerous  successful  examples  of  this 
approach  may  be  found  in  Dekkcr  and  Peters  (1993),  O’Reilly  et  al.  (1998),  Kalio  et  al.  (2001), 
Doxaran  et  al.  (2002),  Brezonik  et  al.  (2005),  Dall’Olmo  and  Gitelson  (2005),  Ouillon  et  al.  (2008), 
Lunctta  et  al.  (2009),  Morel  and  Gentili  (2009).  However,  these  ratios  are  not  free  from  uncertainties 
caused  primarily  by  lighting  geometry,  particulate  scattering  phase  function  p(O),  and  contribution  d\ 
of  the  diffuse  downwelling  irradiancc  to  the  total  (direet+diffuse)  downwelling  irradiance  E&.  To 
overcome  these  uncertainties  in  remote-sensing,  the  use  of  the  Gordon’s  parameter  spectral  ratio  is 
recommended  over  the  application  of  spectral  reflectance  ratios.  In  this  paper,  the  above  inaccuracies 
taking  place  under  using  the  spcetral  reflectance  ratios  of  two  different  types  of  underwater 
reflectances  (irradiance  reflectance  R  and  remote-sensing  reflectance  rn)  are  examined  and  compared 
to  the  G(X)  ratios. 


CONCISE  DESCRIPTION  OF  THE  MAIN  CALCULATIONS 

The  R  and  r n  may  be  defined  in  terms  of  physical  (optical)  quantities  as  follows  (Haltrin,  1998; 
Sokoletsky,  2005;  Sokoletsky  et  al.,  2009a): 


and 


R  =  dFSA  +  (\  —  dF  )PA  , 


(1) 


rn  =—[c/fSA  +  (l-</£)RF],  (2) 

n 

where  SA,  PA,  and  RF  are  the  spherical  albedo,  plane  albedo,  and  reflectance  factor,  respectively. 
Determinations  of  these  quantities  are  given,  for  instance,  in  Hulst  (1974),  Nieodemus  et  al.  (1977), 
Hapke  ( 1 993),  Kokhanovsky  and  Sokoletsky  (2006),  and  briefly  repeated  below.  The  refleetanee  faetor 
is  defined  as  the  ratio  of  the  intensity  of  light  reflected  from  a  given  layer  to  the  intensity  of  light 
reflected  from  the  Lambertian  absolutely  white  surfaee.  The  plane  albedo  is  defined  as  the  integral  of 
the  refleetanee  faetor  RF(//,,  //v ,  (p)\ 


1  zn  i 

PA (//,. )  =  —  {  j  RF(/if ,  p, ,  (p)fjvdpvdcp , 
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(4) 


and  the  spherical  albedo  is  defined  as  the  integral  of  the  plane  albedo: 

i 

SA=2jPA(/i, )//,«///,, 

0 

where  p\  is  the  cosine  of  the  incidence  angle  0,  in  the  medium  (water),  //v  is  the  cosine  of  the  viewing 
angle  <?v  in  the  medium,  and  tp  is  the  azimuthal  angle  between  the  incident  and  scattered  beam 
directions.  Note  that  it  follows  from  the  conservation  energy  law:  SA  =  1  and  also  PA(//j)  =  1  at  any  //j 
for  nonabsorbing  semi-infinite  media  (Sobolev,  1975;  Kokhanovsky,  2004),  although  RF(//j,  pv,  <p) 
may  be  >  1  at  coo  =  1 ,  where  coo  is  the  single-scattering  albedo.  The  d e  may  vary  approximately  from 
1/8  to  1/2  in  the  visible  range  for  different  atmospheric  and  geometrical  conditions  (0,  <  45°)  (Haltrin, 
1997;  Rigollier  et  al.,  2000;  Hojcrslev,  2004). 

Calculations  of  RF  were  performed  following  an  algorithm  and  FORTRAN  code  developed  by 
Mishchenko  et  al.  (1999)  for  the  numerical  solution  of  the  radiative  transfer  equation  (invariant 
imbedding  method).  PA  and  SA  were  derived  by  numerical  integration  following  Eqs.  (3)  and  (4). 
Input  parameters  for  calculation  RF,  PA,  SA,  R  and  r were:  the  given  scattering  phase  function  p(O) 
(with  the  known  baekseattering  probability  B  =  bjb,  where  b  is  the  total  scattering  coefficient),  ruo  and 
Ox.  The  viewing  angle  0 \  was  accepted  as  0°  relative  to  the  zenith  direction  for  all  simulations;  in  this 
particular  ease  a  dependence  of  RF  (and,  hence,  r re)  on  azimuth  ^disappears. 

5  different  p{0)  were  selected  for  simulations  (Kokhanovsky  and  Sokoletsky,  2006;  Sokolctsky  et  al., 
2009b)  so  that  maximally  to  cover  possible  situations  in  natural  waters  from  the  very  clear  oceanic  to 
completely  turbid  inland  waters  (Table  1,  Fig.  1): 


Table  1.  Scattering  phase  functions  used  for  the  study  simulations. 


# 

'ef,  |im 

n 

k 

B 

X 

1 

5 

1.2 

0.01 

0.0087 

0.9583 

2 

0.0163 

0.9377 

3 

2 

1.2 

0.01 

0.0229 

0.9062 

4 

0.175 

1.25 

0.001 

0.0588 

0.6963 

5 

0.116 

1.25 

0.001 

0.1559 

0.5033 

Table  1  shows  input  parameters  used  for  simulation  of  different  p(O):  effective  radius  rc  r  of  the  gamma 
particle  size  distribution  and  refractive  indices  m  =  n-ik\  asymmetry  parameter  g  and  baekseattering 
probability  B  for  the  derived  p{0)  are  also  shown.  Note  that  the  phase  function  #2  is  a  quasi-Foumier- 
Forand-Mobley  p{0)  is  not  directly  connected  with  rc(  and  nr.  this p(O)  has  been  described  separately 
by  Sokoletsky  et  al.  (2009b). 
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Scattering  Angle  (deg) 

Fig.  L  Five  scattering  phase  functions  used  for  the  study  simulations. 


IRRADIANCE  REFLECTANCE  AND  REMOTE-SENSING  REFLECTANCE  SIMULATIONS 


Simulations  of  R  and  r ^  as  a  function  of  G—  0M)BI[\~  tt>o(l-Z?)]  performed  for  5  selected  phase 
functions,  listed  in  the  Table  1 ,  two  incidence  angles  6X  in  the  water  (0°  and  46°)  and  two  values  of  dy 
(1/8  and  1/2)  arc  shown  in  Figs.  2  and  3.  As  it  seems  from  these  figures,  the  curves  may  be  well 
presented  by  the  linear  dependence  in  the  double  logarithmic  scales,  hence,  R  and  r re  may  be  written  as 
power  dependences  on  G : 


R(G,  p(0),  9„  dE)  =  a„ (p(ff),  0„  dE)GMpie'-°'-d'\ 


(5) 


rr,(G,  p(0),  0 „  0V,  <p,  dE)  =  ar(p(6),  6», ,  0V,  <p,  dE)G 


Pr(p(Q),  0;  .  .  <P.(iE  ) 


.(6) 


It  is  important  to  stress  that  irradiancc  reflectance  depends  on  lesser  number  of  parameters  than 
remote-sensing  reflectance.  Further  simplifications  may  be  made  by  applying  the  spectral  reflectance 
ratios: 
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Fig.  2.  Irradiance  reflectance  calculated  for  the  given  p(0)  taken  from  the  Table  I  (the  first  number 
in  parentheses),  0/  (the  second  number,  in  degrees)  and  dy  (the  third  number). 
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Fig.  3.  The  same  as  Fig.  2,  but  for  remote-sensing  reflectance. 
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q  =  rn(G(X ,),  p{0),6;,6v,(p,  dE)  g 
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gw 

GW 


/?r  (  /?( 6? ).  .0V  ,<p.dE) 


=  Q, 


/?r  (/>(<?),  0V .  p,  ) 


(8) 


As  it  can  be  seen  from  Eqs.  (7)  and  (8),  only  Gordon’s  parameter  spectral  ratios  Qg  and  exponents  P r 
and  px  arc  included  in  the  spectral  reflectance  ratio  computations.  The  spectral  differences  for /;(/?),  0t, 
6\,,  and  d\  parameters,  generally  very  small,  were  excluded.  The  values  of  /?r  and  p,  estimated  by  the 
least-squares  method  for  different  p{ 0),  G  <  0.6,  6V  =  0\  and  =1/8  and  d\  =  1/2  arc  shown  as 
functions  of  6\  in  Figs.  4  and  5,  respectively.  The  G  values  used  include  all  known  natural  waters 
(Haltrin  and  Gallegos,  2003).  As  expected,  from  the  optical  behavior  of  RF,  PA  and  SA,  variability  of 
exponents  with  for  R(0,)  is  less  than  for  /*„(#,),  especially  for  small  6\.  Simulations  performed  for  all 
5  p{0)  described  above,  14  angles  G  [0°;  48°],  and  7  values  of  d f  G  [1/8;  1/2],  yielded  values  of  p r 

=  1.04-fc0.02  and  /?r  =  1.05±0.04. 

Variability  of Qr  [Var(fjR)  and  Qr  [Var(i2r)]  with  the  above  parameters  was  estimated  as  the  standard 
deviation  divided  by  the  average  values  of  Qr  and  Qr,  respectively,  versus  Qg  and  d\  using  Eqs.  (7) 
and  (8).  Plots  of  this  variability  (Fig.  6)  show  that  Var(QR)  <  Var(Qr);  the  variability  is  lesser  for 
larger  values  of  d e.  In  any  ease,  the  figure  demonstrates  that  inaccuracies  resulting  from  the  different 
p(0)  and  $  are  rather  small.  It  explains  the  success  of  the  remote-sensing  methods  based  on  spectral 
reflectance  ratios.  This  modeling  has  shown  that  uncertainties  lie  in  the  range  of  [0%;  4.1%]  and  [0%; 

1 1.3%]  for  Qr  and  Qr,  respectively,  if  the  Qg  G  [0.1;  10],  However,  within  the  more  realistic  range  of 
Qg  G  [0.5;  2],  the  same  uncertainties  arc  significantly  smaller:  [0%;  1.3%]  and  [0%;  3.1%]. 

It  is  important  to  stress  that  the  Gordon’s  parameter  spectral  ratios  Qg  (over  of  the  spectral  reflectance 
ratios  considered  in  the  study)  do  not  depend  completely  on  illumination/obscrvation  and  atmospheric 
conditions  and,  hence,  their  using  is  more  preferable  than  QR  and  Qr.  The  spectral  values  of  G(A)  can 
be  directly  derived  from  the  inherent  optical  measurements  of  ocean  optics  sensors  such  as  WetLabs 
ac-9  or  ac-s  spectrophotometers. 
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Fig.  4 .  Irradiance  reflectance  exponent  calculated  for  the  given  p(0)  taken  from  the 
Table  1  (the  first  number  in  parentheses)  and  d e  (the  second  number ). 
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Fig.  5.  The  same  as  Fig.  4 ,  but  for  the  remote-sensing  reflectance  exponent. 
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Fig.  6.  Variations  of  irradiance  reflectance  (solid  lines)  and  remote-sensing  reflectance  ( dotted 
lines)  spectral  ratios  as  a  function  of  the  Gordon 's  parameter  spectral  ratio  and  d  f  (legend). 
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